Multipolar Solutions 

Hernando Quevedo 
Instituto de Ciencias Nucleares, 
Universidad Nacional Autonoma de Mexico 
AP 70543, Mexico, DF 04510, Mexico 
ICRANet, Dipartimento di Fisica, 
Universitd di Roma "La Sapienza" 
1-00185 Roma, Italy 

Abstract 

A class of exact solutions of the Einstein-Maxwell equations is presented which contains infinite 
sets of gravitoelectric, gravitomagnetic and electromagnetic multipole moments. The multipolar 
structure of the solutions indicates that they can be used to describe the exterior gravitational 
field of an arbitrarily rotating mass distribution endowed with an electromagnetic field. The 
presence of gravitational multipoles completely changes the structure of the spacetime because 
of the appearance of naked singularities in a confined spatial region. The possibility of covering 
this region with interior solutions is analyzed in the case of a particular solution with quadrupole 
moment. 
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I. INTRODUCTION 



The problem of describing the gravitational field of astrophysical compact objects is very 
important in general relativity. Indeed, one of the most interesting applications of Einstein's 
general relativity consists in describing the exterior and interior gravitational field generated 
by a specific mass distribution. Accordingly, it must be possible to find exact solutions of 
Einstein's field equations which correctly describe the influence of gravity inside and outside 
the mass distribution. 

This problem is solved in Newtonian gravity by finding the explicit form of the gravita- 
tional potential as a solution of the Poisson and the Laplace equations, respectively. To this 
end, it is usually assumed that the mass distribution satisfies certain symmetry conditions 
that are derived from observations. In particular, one can assume that the mass distribution 
does not depend on time and on the azimuthal angle, i.e. it is stationary and axisymmetric. 
In this case, the expression of the gravitational potential can be written as an infinite series 
each term of which represents a particular multipole moment. 

An additional important aspect is that all known astrophysical compact objects rotate 
with respect to observers situated at infinity where the gravitational field is negligible. 
In Newtonian gravity, however, it is well-known that there is no difference between the 
gravitational potentials of a static and of a rotating (stationary) mass distribution. The 
question arises whether in Einstein's theory of gravity one can find exact solutions that 
take into account the rotation and axial symmetry of the mass distribution. It is therefore 
of importance and interest to describe the relativistic gravitational fields of astrophysical 
compact objects in terms of their multipole moments, in close analogy with the Newtonian 
theory, taking into account their rotation and their internal structure. 

In general relativity, the gravitational field of compact objects must be described by a 
metric ga/s (a, (3 = 0,1, 2, 3) satisfying Einstein's equations, 

Ral3 - -^dapR = ^T^Ta^ , (l) 

in the interior part of the object (Tq,^ ^ 0) as well as outside in empty space {Taj3 = 0), where 
Ta^ represents the energy-momentum tensor of the source of gravity. Since at the classical 
level the most important sources of gravity are the mass distribution and the electromagnetic 
distribution, it is clear that the problem of finding a metric that describes the gravitational 
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field of compact objects can be divided into four related problems. The first one consists 
in finding an exact vacuum solution {Rai3 = 0) that describes the field outside the mass 
distribution. Secondly, the electromagnetic distribution must be considered by solving the 
electrovacuum field equations 



(2) 



where Fap is the Faraday tensor corresponding to a charge distribution Q{x"). As for the 
interior region, it is necessary to propose a model to describe the internal gravitational 
structure of the object. The most popular model is that of a perfect fluid 

Tap = (p + p)UaU(S - QapP , (3) 

with density p(a;"), pressure and 4- velocity m". Finally, the fourth part of the problem 

consists in considering the mass and charge distribution simultaneously, i.e. one needs an 
exact solution to the equations 

Rap - ]^9af}R = Stt T„^(p, p, Q) , (4) 

where a specific model must be proposed for the energy-momentum tensor. Each of the above 
problems is very difficult to solve because Einstein's field equations are highly non-linear. 
Soon after the formulation of Einstein's theory of gravity, the first exterior solution with 

n n 

only a monopole moment was discovered by Schwarzschild In 1917, Weyl [2| showed that 
the problem of finding static axisymmetric vacuum solutions can generically be reduced to 
a single linear differential equation whose general solution can be represented as an infinite 
series. The explicit form of this solution resembles the corresponding solution in Newtonian 
gravity, indicating the possibility of describing the relativistic gravitational field by means 
of multipole moments. In 1918, Lense and Thirring [sl discovered an approximate exterior 
solution which, apart from the mass monopole, contains an additional parameter that can be 
interpreted as representing the angular momentum of the massive body. From this solution 
it became clear that in Einstein's relativistic theory rotation generates a gravitational field 
that leads to the dragging of inertial frames (Lense- Thirring effect). This is the so-called 
gravitomagnetic field which is of especial importance in the case of rapidly rotating compact 
objects. The case of a static axisymmetric solution with monopole and quadrupole moment 



was analyzed in 1959 by Erez and Rosen 



J] by using spheroidal coordinates which are 



specially adapted to describe the gravitational field of non-spherically symmetric bodies. 
The exact exterior solution which considers arbitrary values for the angular momentum was 
found by Kerr only in 1963. The problem of finding exact solutions changed dramatically 
after Ernst [6] discovered in 1968 a new representation of the field equations for stationary 
axisymmetric vacuum solutions. In fact, this new representation was the starting point to 
investigate the Lie symmetries of the field equations. Today, it is known that for this special 
case the field equations are completely integrable and solutions can be obtained by using 
the modern solution generating techniques [7| . In thi^ work, we will analyze a particular 
class of solutions, derived by Quevedo and Mashhoon [8] in 1991, which in the most general 
case contains infinite sets of gravitational and electromagnetic multipole moments. 

As for the interior gravitational field of compact objects, the situation is more compli- 
cated. There exists in the literature a reasonable number of interior spherically symmetric 
solutions {9 1 which can be matched with the exterior Schwarzschild metric. Nevertheless, 
a major problem of classical general relativity consists in finding a physically reasonable 
interior solution for the exterior Kerr metric. Although it is possible to match numerically 
the Kerr solution with the interior field of an infinitely tiny rotating disk of dust 10|, such 
a hypothetical system does not seem to be of relevance to describe astrophysical compact 
objects. It is now widely believed that the Kerr solution is not appropriate to describe 
the exterior field of rapidly rotating compact objects. Indeed, the Kerr metric takes into 
account the total mass and the angular momentum of the body. However, the moment of 
inertia is an additional characteristic of any realistic body which should be considered in 
order to correctly describe the gravitational field. As a consequence, the multipole moments 
of the field created by a rapidly rotating compact object are different from the multipole 



moments of the Kerr metric 



ll| . For this reason a solution with arbitrary sets of multipole 



moments, such as the one presented in this work, can be used to describe the exterior field 
of arbitrarily rotating mass distributions. 



II. LINE ELEMENT AND FIELD EQUATIONS 

Although there exist in the literature many suitable coordinate systems, stationary ax- 
isymmetric gravitational fields are usually described in cylindric coordinates {t, p, z, ip). Sta- 
tionarity implies that t can be chosen as the time coordinate and the metric does not depend 



4 



on time, i.e. dgap/dt = 0. Consequently, the corresponding timelike Killing vector has the 
components 5f . A second Killing vector field is associated to the axial symmetry with re- 
spect to the axis p = 0. Then, choosing (f as the azimuthal angle, the metric satisfies the 
conditions dgajs/dLp = 0, and the components of the corresponding spacelike Killing vector 
are 6^. 

Using further the properties of stationarity and axial symmetry, together with the vacuum 
field equations, for a general metric of the form Qa/s = Qapip, z), it is possible to show that 
the most general line element for this type of gravitational fields can be written in the 
Weyl-Lewis-Papapetrou form as j^, I2I ] 
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ds^ = f{dt~ujd^f - f-^[e^^{dp'' + dz'')+ p'd^^] , (5) 

where /, u and k are functions of p and z, only. After some rearrangements which include 
the introduction of a new function Vt = Vt{p, z) by means of 

pdpVL = fd;,u , pd^n = -fdpuj , (6) 

the vacuum field equations Rap = can be shown to be equivalent to the following set of 
partial differential equations 

Upipdpf) + dlf + jiidp^f + id.^f - idpff - [djf] = , (7) 
-^dp{pdp9) + din - J {dpf dpQ + dj d,n) = , (8) 

= {f2 [(^ff + (^p^)' - (d-Jf - {d^m ' (9) 
d,k = ^ {dpf dj + dpn d-M) . (10) 

It is clear that the field equations for k can be integrated by quadratures, once / and Vt are 
known. For this reason, the equations ([7]) and ([8]) for / and are usually considered as the 
main field equations for stationary axisymmetric vacuum gravitational fields. 

It is interesting to mention that the main field equations can be obtained from a La- 
grangian in the following way. The Einstein-Hilbert Lagrangian Ceh = V~9^ ^i^® 
element (jS]) with the auxiliary function Q, as defined in Eq.(l6]), can be written as 

Ceh = ^ [{dpfY + {djf + {dp^Y + [d^^lf] , (11) 
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where the terms containing second order derivatives have been ehminated by neglecting the 
total divergence terms, and a Legendre transformation has been applied for the "cyclic" 
functions Q and k 15[. Then, the variation of this Lagrangian density with respect to / and 



Q generates the main field equations ([7j) and ([H]). 
A. Generalized harmonic maps 

An alternative differential geometric interpretation of stationary axisymmetric gravita- 
tional fields can be explored by using the concepts of harmonic maps [l^ as follows. 

Consider two (pseudo-)Riemannian manifolds (M, 7) and {N,G) of dimension m and n, 
respectively. Let x° and be coordinates on M and A^, respectively. This coordinatization 
implies that in general the metrics 7 and G become functions of the corresponding coordi- 
nates. Let us assume that not only 7 but also G can explicitly depend on the coordinates 
x", i.e. let 7 = 7(x) and G = G{X,x). A smooth map X : M ^ N will be called an 
(m — )• n) —generalized harmonic map if it satisfies the Euler-Lagrange equations 

-^d, (^/Wb^'daX'^) + d^X-d^X^ + ^^^7^^' daX^ d^Gx, = , (12) 

which follow from the variation of the generalized action 

S = j d^x^/W\r\x)daX>^dhX''G^,{X,x) , (13) 

with respect to the fields X^. Here the Christoffel symbols, determined by the metric G^^, 
are calculated in the standard manner, without considering the explicit dependence on x. 
Notice that the presence of the term G^i^{X, x) in the Lagrangian density takes into account 
the "interaction" between the base space M and the target space N. This interaction leads 
to an extra term in the motion equations, as can be seen in ( fT2l) . which is important in 
order to recover the correct field equations in the case of gravitational fields. Moreover, 
this interaction affects the conservation laws of the physical systems we attempt to describe 
by means of generalized harmonic maps. To see this explicitly we calculate the covariant 
derivative of the generalized Lagrangian density 

^ = VH7"''(a;)9aX^9,X^G^,(X,x) , (14) 

and replace the result in the corresponding motion equations (|T2|) . Then, the final result 



can be written as 



V.T/ + i|^ = 0. (15) 



where Tab represents the canonical energy-momentum tensor 

T., = ^ , = ^G,. (^I'^'daX^ d,X^ - Ui^^'d^X^ daX'^'^ . (16) 

The standard conservation law (Vf,T^^ = 0) is recovered only when the Lagrangian does 
not depend explicitly on the coordinates of the base space. Even if we choose a flat base 
space = r/ab, the explicit dependence of the metric of the target space G^iy(^X, x) on 
X generates a term that violates the standard conservation law. This term is due to the 
interaction between the base space and the target space which, consequently, is one of the 
main characteristics of the generalized harmonic maps. 

Consider a (2 — )■ 2)— generalized harmonic map. Let = {p,z) be the coordinates on 
the base space M, and = (/, Q) the coordinates on the target space A^. In the base 
space we choose a flat metric and in the target space a conformally fiat metric, i.e. 

lab = Sab and Gf,^ = ^6^,^ (a, 6 = 1, 2; fi,u = 1,2). (17) 

A straightforward computation shows that the generalized Lagrangian (fT4|) coincides with 
the Lagrangian ( ITTl) for stationary axisymetric fields, and that the equations of motion ( !T2|l 
generate the main field equations (JTj) and (jH]). Moreover, if we calculate the components of 
the energy- momentum tensor Tab = SC/6'~f"'^, we obtain 

Tpp = -T.. = ^ [{d.fY + {d,nr - (djr - [d.nr] , (is) 

Tp. = ^ {d,f dj + dpQ d^Q) . (19) 

This tensor is traceless due to the fact that the base space is 2-dimensional. It satisfies the 
generalized conservation law f lT5|) on-shell: 

^ - ^ = . (21) 

dp dz 

Incidentally, the last equation coincides with the integrability condition for the metric func- 
tion k, which is identically satisfied by virtue of the main field equations. In fact, as can be 
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seen from Eqs.( l9|T0|) and ( !T8l[T9i) . the components of the energy-momentum tensor satisfy 
the relationships Tpp = dpk and Tpz = dzk, so that the conservation law fl2Tl) becomes an 
identity. Although we have eliminated from the starting Lagrangian f llip the variable k by 
applying a Legendre transformation on the Einstein-Hilbert Lagrangian (see [isl for details) 
for this type of gravitational fields, the formalism of generalized harmonic maps seems to 
retain the information about k at the level of the generalized conservation law. 

The above results show that stationary axisymmetric spacetimes can be represented as 
a (2 —7- 2)— generalized harmonic map with metrics given as in (1171) . It is also possible to 
interpret the generalized harmonic map given above as a generalized string model. Although 
the metric of the base space M is Euclidean, we can apply a Wick rotation t = ip to obtain 
a Minkowski-like structure on M. Then, M represents the world-sheet of a bosonic string 
in which r is measures the time and z is the parameter along the string. The string is 
"embedded" in the target space N whose metric is conformally flat and explicitly depends 



on the time parameter r. For more details see 17l |. 



III. THE STATIC SOLUTION 

Let us consider the special case of static axisymmetric fields. This corresponds to metrics 
which, apart from being axially symmetric and independent of the time coordinate, are 
invariant with respect to the transformation (f — )■ —(f (i.e. rotations with respect to the axis 
of symmetry are not allowed). Consequently, the corresponding line element is given by 

ds^ = fdt^-f-^[e^''{dp^ + dz^)+p^d^^] , (22) 

and the field equations can be written as 

dl^P + -dp^ + a^V- = , / = exp(2^) , (23) 

dpk = p [{dp^)^ - id^^'f] , d,k = 2pdp^ <9,V . (24) 

We see that the main field equation f l23|) corresponds to the linear Laplace equation for the 
metric function ib. The general solution of Laplace's equation is known and, if we demand 
additionally asymptotic flatness, we obtain the Weyl solution which can be written as [2, Il4j | 

oo 

^ = E , , „-^. Pn(cOS 9) , COS ^ = -^^== , (25) 



where a„ (n = 0, 1, ...) are arbitrary constants, and P„(cos^) represents the Legendre polyno- 
mials of degree n. The expression for the metric function 7 can be calculated by quadratures 
by using the set of first order differential equations fl24l) . Then 

^^^^ (n + m + 2)(p2 + ^2)^ 

Since this is the most general static, axisymmetric, asymptotically flat vacuum solution, it 
must contain all known solutions of this class. In particular, one of the most interesting 
special solutions which is Schwarzschild's spherically symmetric black hole spacetime must 
be contained in this class. To see this, we must choose the constants a„ in such a way that 
the inflnite sum (1251) converges to the Schwarzschild solution in cylindric coordinates. But, 
of course, this representation is not the most appropriate to analyze the interesting physical 
properties of Schwarzschild's metric. 

In fact, it turns out that to investigate the properties of solutions with multipole moments 
it is more convenient to use prolate spheroidal coordinates (t, x, y?) in which the line 
element can be written as 



ds"" = fde - ^ 



2 



e'\x' - y') ( ^ + i)(i - y^) V 



(27) 



where 



x='-^^. (x^>i), y='-hr-^ (^'^1) (28) 



= + (;z ± a) , a = const , (29) 

and the metric functions /, and k depend on x and y, only. In this coordinate system, the 
main fleld equation becomes 

[{x^ - + [(1 - y^)ijy]y = , / = exp(2V^) , (30) 

and the general static solution which is also asymptotically flat can be expressed as 

00 

ip = J](-l)"+^g„P„(y)Q„(x) , g„ = const (31) 

n=0 

where Pn{y) are the Legendre polynomials, and Qn{x) are the Legendre functions of second 
kind. In particular, 

Po = i, Pi = y, P2 = ^(3i/'-i) 
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1 x+1 1 x+1 
Qq = - In , Q^ = -X In r - 1 



2 x-1 



Q2 = l(3a;2_i)ij^^_3 
2^ ' x~l 2 



2 x-1 
3 

— -a; 



The corresponding function k can be calculated by quadratures and its general expression 
has been explicitly derived in jisl . 

The most important special cases contained in this general solution are the Schwarzschild 
metric 

ij = -qoPo{y)Qo{x) = - In , k = -In . (32) 

Indeed, the coordinate transformation 

r — m 



y = cos 9 , X 



a 



(33) 



transforms the line element (1271) into 



ds^ = fdt^ - i 



-,2k 



2m I m^—(T^ cos^( 



2m I m^— (7^ 



f 1 - — + ) r2 sin^ 



l-iHi + i 

r 

2 



(34) 



and the metric functions (1321) with go = 1 into 



2 ~ r / ' 2 \ 1 2m I — (T^ cos^ ' 



(35) 

which, when inserted in the above line element with a = m, lead to the standard 
Schwarzschild solution in spherical coordinates 



1 ]df 



2m 



r\de' + sin' 



r J 1 
fl 

Moreover, the Erez-Rosen metric is obtained in the special case 



(36) 



iP = -P,{y)Qo{x) - q2P2{y)Q2{x) 

= iinf^^ +ig,(3y2_i) 
2 \x+lj 2^^^ ^ ' 



i(3x2-l)ln 



and 



k= ^{1 + q,)Hn ^ ^ 



— y'^ 



l<i2{i-y') 



x\n ( ^ 1 1 +2 

x + 1 



+ ^^2^(1-1/^) 



3 



(37) 



x^ + 4y'' - Ox'^y^ \-x [x' + 7y- 9x'y' In 



X — 1 
x + 1 



+ ^(X^ - 1)(X2 + y2 _ 9a;2^2 _ ^) 1^2 



(38) 
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In the last case, the constant parameter q2 turns out to determine the quadrupole moment. 
In general, the constants g„ represent an infinite set of parameters that determines an infinite 
set of mass multipole moments. In fact, using the Geroch-Hansen 19|, |20| definition, one 
can prove that the relativistic multipole moments can be expressed as 

11} 

Mn = Nn + Rn, Nn= ('^ )" + 1 ) ! ! ^""^ ' ^ = 0' 1' 2' (39) 

where Nn are the Newtonian multipole moments which have been calculated by using the 
coordinate invariant method proposed by Ehlers in j21|. Moreover, the second term i?„ 
represents the relativistic corrections 

i?o = = i?2 = , i?3 = -^fT^iVi , i?4 = ~(y^N2 - ^-aNl , .... (40) 

which in general can be determined in terms of lower Newtonian moments, i. e., Rn = 

RniNn-2,Nn-3,-,No). 

The metric function given in Eq.f pT]) . together with the corresponding function k derived 
in 18|, represents the most general static axisymmetric solution of vacuum Einstein's equa- 
tions in prolate spheroidal coordinates. Its multipolar structure represented by the infinite 
set of parameters g„ can be used to describe the exterior gravitational field of any static 
mass distribution that preserves the axial symmetry. 

The simplest metric contained in this class of solutions corresponds to the Schwarzschild 
spacetime which describes the exterior gravitational field of a black hole of mass m. Ac- 
cording to the black hole uniqueness theorems the Schwarzschild spacetime represents the 
only static black hole, i. e., it is the only solution with a singularity covered by an event 
horizon. It then follows that all the multipolar solutions with multipoles higher than the 
monopole one must be characterized by the presence of naked singularities. This has been 
shown explicitly for the Erez-Rosen solution in 22[ by using a numerical approach due to 
the complexity of the resulting curvature invariants. To be able to perform an analytical 
investigation it is necessary to consider a simpler metric. To this end, we use the following 
property of the field equations for static fields. If ipQ and fco represent an exact static solution 
of the field equations fl23|) and (^^, then the functions 6ipo and S^ko, with 6 = const are also 
solutions to the same field equations. This property was first discovered by Zipoy [23[ and 



Voorhees 



24| . Consider then the solution 

(5 , X — 1 , (5^ , 

- In , k = — In — 

2 x + 1 ' 2 x^-y 



= - In — ^ , k = — In _ ^ , (41) 
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which represents the simplest generahzation of the Schwarzschild metric ( 132|) with go = 1- 
Then, introducing spherical coordinates by means of the relations f l33|) . and choosing the 
parameters a = m and 5 = 1 — g, the resulting solution can be written as 



r y 

r 



This solution is axially symmetric and reduces to the spherically symmetric Schwarzschild 
metric in the limit g — )■ 0. It is asymptotically flat for any finite values of the parameters 
m and g. To find the physical meaning of these parameters, we calculate the multipole 



moments of the solution by using the invariant definition proposed by Geroch 19|. The 
lowest mass multipole moments M„, n = 0,1, . . . are given by 

Mo = (l-g)m, M2 = — g(l-g)(2-g) , (43) 

whereas higher moments are proportional to mq and can be completely rewritten in terms 
of Mq and M2. This means that the arbitrary parameters m and g determine the mass 
and quadrupole which are the only independent multipole moments of the solution. In the 
limiting case g = only the monopole Mq = m survives, as in the Schwarzschild spacetime. 
In the limit m = 0, with g 7^ 0, all moments vanish identically, implying that no mass 
distribution is present and the spacetime must be flat. This can be seen also at the level of 
the curvature which vanishes in the limiting case m — )■ 0. This means that, independently of 
the value of g, there exists a coordinate transformation that transforms the resulting metric 
into the Minkowski solution. From a physical point of view this is an important property 
because it means that the parameter g is related to a genuine mass distribution, i.e., there is 
no quadrupole moment without mass. Furthermore, notice that all odd multipole moments 
are zero because the solution possesses an additional reflection symmetry with respect to 
the equatorial plane. 

We conclude that the above metric describes the exterior gravitational field of a static 
deformed mass. The deformation is described by the quadrupole moment M2 which is 
positive for a prolate mass distribution and negative for an oblate one. Notice that in order 
to avoid the appearance of a negative total mass Mq the condition g < 1 must be satisfied . 

To investigate the structure of possible curvature singularities, we consider the 
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Kretschmann scalar K = Raises R"^^^ ■ A straightforward computation leads to 
_ 16m2(l - g)2 (r^ - 2mr + sin^ ^) V-4g-i 



with 

L(r, ^) = 3(r — 2m + qmY{r'^ — 2mr + sin^ 9) 

— g(2 — g)m^ sin^ 9[q{2 — q)m'^ + 3(r — m)(r — 2m + qm)] . (45) 

In the limiting case g = 0, we obtain the Schwarzschild value K = A8m^/r^ with the only 
singularity situated at the origin of coordinates r — )■ 0. In general, one can show that 
the singularity at the origin, r = 0, is present for any values of q. Moreover, an additional 
singularity appears at the radius r = 2m which, according to the metric ( H3|) . is also a horizon 
in the sense that the norm of the timelike Killing tensor vanishes at that radius. Outside the 
hypersurface r = 2m no additional horizon exists, indicating that the singularities situated 
at the origin and at r = 2m are naked. Moreover, for values of the quadrupole parameter 
within the interval 

g G (l - v^, 1 + \{0} (46) 

a singular hypersurface appears at a distance 

r± = m(l ± cos6') (47) 

from the origin of coordinates. This type of singularity is always contained within the naked 
singularity situated at the radius r = 2m, and is related to a negative total mass Mq for 
g > 1. Nevertheless, in the interval g G (1 — a/3/2, 1]\{0} the singularity is generated 
by a more realistic source with positive mass. This configuration of naked singularities is 
schematically illustrated in Fig. [H 

Another important aspect related to the presence of naked singularities in multipolar 
solutions is the problem of repulsive gravity. In fact, it now seems to be established that 



naked singularities can appear as the result of a realistic gravitational collapse [3l| and that 
naked singularities can generate repulsive gravity. Currently, there is no invariant definition 
of repulsive gravity in the context of general relativity, although some attempts have been 



made by using invariant quantities constructed with the curvature of spacetime 25l-l27|. 



Nevertheless, it is possible to consider an intuitive approach by using the fact that the 
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(b) 



FIG. 1: Structure of naked singularities in a spacetime with quadrupole parameter q. For any 
value q ^ 0, there exists at least two naked singularities situated at r = and r = 2m as shown 
in plot (b) with solid curves. Furthermore, if the quadrupole parameter is contained within the 



plot (c). The limiting case of a Schwarzschild spacetime {q = 0) with a singularity at the origin of 
coordinates surrounded by the horizon (dashed curve) situated at r = 2m is illustrated in plot (a) . 

motion of test particles in static axisymmetric gravitational fields reduces to the motion in 
an effective potential. This is a consequence of the fact that the geodesic equations possess 
two first integrals associated with stationarity and axial symmetry. The explicit form of the 
effective potential depends also on the type of motion under consideration. In the case of a 
massive test particle moving along a geodesic contained in the equatorial plane {9 = vr/2) 
of the Zipoy-Voorhees spacetime ( H3l) . one can show that the effective potential reduces to 



where L is constant associated to the angular momentum of the test particle as measured by 
a static observer at rest at infinity. This expression shows that the behavior of the effective 
potential strongly depends on the value of the quadrupole parameter q. This behavior is 
illustrated in Fig. |2l 

Whereas the effective potential of a black hole corresponds to the typical potential of an 
attractive field, the effective potential of a naked singularity is characterized by the presence 
of a barrier which acts on test particles as a source of repulsive gravity. We can show that 
in general in a mass distribution the presence of static multipoles higher than the monopole 
one leads to the appearance of naked singularities in which effects associated with repulsive 
gravity can be found. It is therefore clear that the mass quadrupole and higher multipoles 
can be considered as sources of naked singularities in general relativity. 



interval q £ 



1 — a/3/2, 1 + \/3/2 ] \{0}, two additional naked singularities appear as depicted in 




(48) 
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FIG. 2: The effective potential for tlie motion of timelike particles. Plot (a) shows the typical 
behavior of the effective potential of a black hole configuration with q = 0. The case of a naked 
singularity with q = 1/2 is depicted in plot (b). 

IV. STATIONARY SOLUTIONS 

The solution generating techniques can be applied, in particular, to any static seed 

solution in order to obtain the corresponding stationary generalization. One of the most 
powerful techniques is the inverse scattering method developed by Belinski and Zakharov 
. In this section we use a particular case of the inverse scattering method which is known 



as the Hoenselaers-Kinnersley-Xanthopoulos (HKX) transformation to derive the stationary 
generalization of the general static solution in prolate spheroidal coordinates presented in 
the last section. 



A. Ernst representation 



In the general stationary case (u; 7^ 0), the line element in prolate spheroidal coordinates 
is given as 



ds'^ = f{dt — udipY 



a 



f 



dx^ dy 



e-(x^ - y^) ( ^ + ^ ) + (x^ - 1)(1 - V 



(49) 
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where all the metric functions depend on x and y only. It turns out to be useful to introduce 
the the complex Ernst potentials 

E = f + in, e = ^^, (50) 

where the function Vl is now determined by the equations 

a{x^ - = puy , a(l - y^)ny = -fu^ . (51) 

Then, it is easy to show that the main field equations can be represented in a compact and 
symmetric form as 

(er - 1) {[{x' - + [(1 - y')Qy] = 2ei(x^ - + (1 - y')Q , (52) 

where the asterisk represents complex conjugation. Notice that in the case of static fields the 
Ernst potentials become real and the above equation generates a linear differential equation 
for %l) = (1/2) In/. It is easy to see that the equation ( 152|) is invariant with respect to the 
transformation x -H- y. Then, since the particular solution 

1 1 - 1 

^ = -^Q = 0-^uj = 0^k=-\n — (53) 

X 2 x^ — 

represents the Schwarzschild spacetime, the choice = 2/ is also an exact solution. Fur- 
thermore, if we take the linear combination C,~^ = cix + C2y and introduce it into the field 
equation (l30i) . we obtain the new solution 



^ ^ = —X + i—y , cr = \/m? — , (54) 
m m 



which corresponds to the Kerr metric in prolate spheroidal coordinates. The corresponding 
metric functions are 

cV + d^y^-l _ (ca; + 1)(1 -|/^) 

/ — 7 , -1X9 I — vT^ ■> w — 2a 



cx + 1)^ + (Py"^ ' c^x^ + (Py"^ — 1 ' 

1 / cV + dV- r 

2 V C2(x2-y2 



^ = oln( — ^-^--^ ) , (55) 



where 



c=-, rf=-, ^ + d^ = l. (56) 

m m 



In the case of the Einstein-Maxwell theory, the main field equations can be expressed as 

(a* - J" J"* - = 2(rV^ - rVJ')Vi , (57) 

16 



- - 1) V'J^ = 2(f - rVJ=)VT (58) 

where V represents the gradient operator in prolate spheroidal coordinates. Moreover, the 
gravitational potential ^ and the electromagnetic T Ernst potential are defined as 

^ 1 + / + ' 1 + / + ^ ' 

The potential $ can be shown to be determined uniquely by the electromagnetic potentials 
Ai and Furthermore, one can show that if is a vacuum solution, then the new potential 



i = CoVT^ (60) 

represents a solution of the Einstein-Maxwell equations with effective electric charge e. This 
transformation is known in the literature as the Harrison transformation 29|]. Accordingly, 
the Kerr-Newman solution in this representation acquires the simple form 



e= /Vf , e = ^, a=y/m^-a^-QK (61) 
—X + t—v m 

In this way, it is very easy to generalize any vacuum solution to include the case of electric 
charge. More general transformations of this type can be used in order to generate solutions 



with any desired set of gravitational and electromagnetic multipole moments 30|. 



V. THE GENERAL SOLUTION 

If we take as seed metric the general static solution in prolate spheroidal coordinates with 
an arbitrary Zipoy-Voorhees parameter, 

oo 

^ = <5^(-l)"g„P„(y)Q„(x) , (62) 

n=0 

the application of two HXK transformations generates a stationary solution with an infi- 
nite number of gravitoelectric and gravitomagnetic multipole moments. The HKX method 
generates a new Ernst potential ^ which can be written as 

e= ^"^ + ^^^^^''' + "~+^^- (l-eg + ,g)V^ $ = ^^L±l^, (63) 



where 



^ = Y.{-irqnPn{y)Qn{x) (64) 



n=l 
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with 



A 

fi -- 



a± = {x± - ± (1 + A/i)] , 

6± = (x±l)^-i[y(A + /i)T(A-/i)] , 



(65) 
(66) 

(67) 
(68) 



and 



n± 



fill 



+ i^«(l/)Qn-l(x) 



n-1 



5^(±l)'^P„_fe(l/) [g„_fc+i(x) - g„-fc-i(a;)] . 



(69) 



k=l 



Here Pn(y) and Qn{x) represent the Legendre polynomials and functions of second kind, 
respectively. The constant parameters Cq, go, ai, 0^2, Qn, and 5 determine the gravitational 
and electromagnetic multipole moments. The metric functions / and u can be obtained 
from the definitions of the Ernst potentials whereas the function k can be calculated by 
quadratures once / and u are known. 

In general, this solution is asymptotically fiat and free of singularities along the axis 
of symmetry, y = 1, outside certain region situated close to the origin of coordinates. 
The sets of infinite multipole moments can be chosen in such a way as to reproduce the 
shape of ordinary axially symmetric compact objects. One of the most interesting solutions 
contained in this family is the one with non-vanishing parameters go = 1)^2 = oi\ = 
0^2 = (c — m)/a, where m and a are new constants. In this case, the solution possesses the 
following independent parameters: m, a, 6, and q. In the limiting case a = 0, a = 0, g = 
and 6 = 1, the only independent parameter is m and the Ernst potential ( l63l) determines the 
Schwarzschild spacetime. Moreover, for a = a = and g = we obtain the Ernst potential 



of the Zipoy-Voorhees (ZV) [23|, |2J] static solution which is characterized by the parameters 
m and S. Furthermore, for a = a = and 6 = 1, the resulting solution coincides with the 
Erez- Rosen (ER) static spacetime The Kerr metric is also contained as a special case 
for q = and 6 = 1. The physical significance of the parameters entering this particular 
solution can be established in an invariant manner by calculating the relativistic Geroch- 



Hansen 



20| multipole moments. We use here the procedure formulated in Ref. 
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which allows us to derive the gravitoelectric M„ as well as the gravitomagnetic J„ multipole 
moments. A lengthy but straightforward calculation yields 

M2fe+i = J2fe = 0, A; = 0,1, 2,... (70) 

Mo^m + a{S- 1) (71) 
M2 = ^a^Sq - \a%S^ - 35^ - 45 + 6) - ma'^S{S - 2) - 3mV(5 - 1) - m\ (72) 

Ji = ma + 2a(7((5 - 1) , (73) 

4 

2 



J3 = ^aa^Sq 
15 



—a 



-a^{5^ -35^-5 + 3)+ ma^{3S^ -66 + 2) + 4mV(5 - 1) + 
3 



(74) 



The even gravitomagnetic and the odd gravitoelectric multipoles vanish identically because 
the solution possesses and additional reflection symmetry with respect to the hyperplane 
y — which can be interpreted as the equatorial plane. Higher odd gravitomagnetic and even 
gravitoelectric multipoles can be shown to be linearly dependent since they are completely 
determined in terms of the parameters m, a, q and 5. Prom the above expressions we 
see that the ZV parameter 5 enters exphcitly the value of the total mass Mq as well as 
the angular momentum Ji of the source. The mass quadrupole M2 can be interpreted 
as a nonlinear superposition of the quadrupoles corresponding to the ZV, ER and Kerr 
spacetimes. A generalization of the Kerr metric which includes an arbitrary quadrupole 
moment is obtained by imposing the condition 5—1. The resulting multipoles are 

M2fe+i = J2fe = 0, A; = 0,1, 2,... (75) 
2 / an 

Mo = m , M2 = -ma^ + —qm^ 1 -) , ... (76) 

15 \ J 

4 f a'^\ 

Ji = ma , J3 = -ma^ + —qm^a 1 ] , .... (77) 

15 \ J 

It is interesting to note that this particular exact solution in the limit a ^ m leads to the 
spacetime of an extreme Kerr black hole, regardless of the value of the quadrupole parameter 

In all the above special solutions the electromagnetic field vanishes identically. One can 
easily obtain the corresponding electrovacuum generalizations by assuming that eo 7^ and 
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qq 7^ 0. The computation of the respective electromagnetic multipole moments can be 
performed in an invariant manner and the result can be expressed as 



En = eoM„ , Hn = goJn ■ (78) 

This means that the charge distribution resembles the mass distribution. The electric mo- 
ments vanish identically if no mass distribution exists. This result is in accordance with our 
physical intuitive interpretation of a charge distribution. The magnetic moments turn out 
to be proportional to the gravitomagnetic multipoles, with no magnetic monopole. This is 
a physical reasonable result in the sense that the magnetic field is generated by the motion 
of the charge distribution, i. e., in the present case, by the rotation of the compact object. 

It is worth noticing that in all the above special solutions we assumed that ai = a2 
and obtained generalizations of the Kerr metric with arbitrary quadrupole moment. More 
general solutions can be obtained by relaxing this condition. Consider, for instance, the 
special solution with 6 = 1, qo = 1, qi = for i = 1, 2, and 



— — r > «2 = r ' ^ = — 2 " ' ^ = 



The resulting potential corresponds to the charged Kerr-Taub-NUT spacetime with total 
charge Qo = meo, where I is the Taub-NUT parameter. 

The analysis of the general stationary solution turns out to be very complicated because 
of its mathematical complexity. The special case of the Kerr metric with only an additional 



quadrupole parameter q2 was recently analyzed in detail in j32j. It was shown that the 
presence of the quadrupole parameter completely changes the structure of spacetime due 
especially to the fact that a naked singularity appears that affects the geometric structure 
of the ergosphere and the motion of test particles around the mass distribution. We can 
expect that similar effects will occur if higher multipole moments are taken into account. 



A. An interior solution 



A major problem in general relativity is to find physically meaningful solutions that can 
be matched with exterior exact solutions. In the context of multipolar solutions, one can 
say that only a few interior Schwarzschild solutions are known which can be considered as 
physically reasonable. The search for an interior solution that could be matched with the 
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Kerr metric is still an open problem. Only recently it was proposed to use the quadrupole 
moment as a parameter that introduces an additional degree of freedom into the differential 
equations which determine the internal structure of the mass distribution 33| . To illustrate 
the method we consider the simplest generalization of the Schwarzschild spacetime with 
quadrupole moment given in Eq.f l43p . In the search for the corresponding interior solution, 
we found that an appropriate form of the line element can be written as 



where 

g2fco = (^2 _ + cos^ e)e^^^''''^^ , (81) 

and / = f{r,6), h = h{r), and fi = fi{r,6). This line element preserves axial symmetry 
and staticity. In general, in order to solve Einstein's equations with a perfect fluid source, 
the pressure and the energy must be functions of the coordinates r and 6. However, if we 
assume that p = const, the complexity of the corresponding differential equations reduces 
drastically: 

Pr = ~{p + p)j, Pe = ~{p + p)j, (82) 
A^rr = i^Pee + KfJ'T - 327rp— , (83) 

Moreover, the function k turns out to be determined by a set of two partial differential equa- 
tions which can be integrated by quadratures once / and /i are known. The integrability 
condition of these partial differential equations turns out to be satisfied identically by virtue 
of the remaining field equations. It is then possible to perform a numerical integration by 
imposing appropriate initial conditions. In particular, we demand that the metric functions 
and the pressure are finite at the axis. It turns out to be possible to find numerical solu- 
tions for the metric functions and the thermodynamic variables. In particular, the pressure 
behaves as shown in Figj3l 

It can be seen that the pressure is finite in the entire interior domain, and tends to zero 
at certain hypersurface R{r, 9) which depends on the initial value of the pressure on the axis. 
Incidentally, it turns out that by increasing the value of the pressure on the axis, the "radius 
function" i?(r, 9) can be reduced. Furthermore, if we demand that the hypersurface /2(r, 9) 
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FIG. 3: Plot of the inner pressure as a function of the spatial coordinates. 

coincides with the origin of coordinates, the value of the pressure at that point diverges. 
From a physical point of view, this is exactly the behavior that is expected from a physically 
meaningful pressure function. 

This solution can be used to calculate numerically the corresponding Riemann tensor 
and its eigenvalues. As a result we obtain that the solution is free of singularities in the 
entire region contained within the radius function R{r,6). To perform the matching with 



the exterior metric (143|) we applied the method formulated in 
properties of the eigenvalues of the curvature tensor. 



which uses the invariant 



VI. CONCLUSIONS 



We presented in this work a class of electrovacuum solutions of Einstein-Maxwell equa- 
tions which can be used to describe the exterior gravitational field of a rotating distribution 
of mass endowed with an electromagnetic field. This class of solutions is characterized by 
different sets of arbitrary parameters which determine the multipolar structure of the grav- 
itational source. An important consequence of the presence of multipoles higher than the 
monopole is that the structure of the spacetime completely changes due to the appearance 
of naked singularities. In all the cases we investigated, the curvature singularities are sit- 
uated inside the horizon which becomes the outermost naked singularity. This means that 
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in principle it should be possible to find an interior solution that covers the entire spatial 
region where the naked singularities exist. In particular, we found numerically an inner 
solution that can be matched with an exact exterior metric with a particular quadrupole 
moment. The entire spacetime is shown to be free of singularities so that the entire man- 
ifold is well-defined in terms of solutions of Einstein's equations. However, the problem of 
finding interior solutions taking into account the rotation of the gravitational source remains 
an open problem. We believe that the multipolar structure of the solutions presented here 
could be used as additional degrees of freedom to search for more realistic inner solutions. 
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